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Abstract. We study the generalization of the Biot-Savart law from electrodynamics in 
the presence of curvature. We define the integral operator BS acting on all vector fields 
on subdomains of the three-dimensional sphere. By doing so, we establish a geometric 
setting for electrodynamics in positive curvature. When applied to a vector field, the Biot- 
Savart operator behaves like a magnetic field; we demonstrate that Maxwell's equations 
hold for this context. In particular, for vector fields that act like currents, the curl operator 
is a left inverse to BS; thus the Biot-Savart operator is important in the study of curl 
eigenvalue energy-minimization problems in geometry and physics. We show that the Biot- 
Savart operator is self-adjoint and bounded. In all instances, the formulas we give are 
geometrically meaningful: they are preserved by orientation-preserving isometries of the 
three-sphere. 



1. Introduction 

The Biot-Savart law in electrodynamics calculates the magnetic field B arising from a 
current flow V in a smoothly bounded region Vt of M 3 as 

(1) B(V)(y) = ^ [ V(x) x dx. 

^ Jn \y-x\ 6 

Taking the curl of B recovers the flow V, provided there is no time-dependence for this 
system. The Biot-Savart law can be extended to an operator BS which acts on all smooth 
vector fields V defined in fi. Cantarella, DeTurck, and Gluck studied the properties of this 
operator and its applications in J6]|. 

The Biot-Savart operator is closely connected to the helicity of a vector field, which 
measures the average linking of its flowlines. Helicity can be computed as the L 2 inner 
product (V,BS(V)). Woltjer [|20]| first applied helicity to astrophysics in 1958. Moffatt 
[fT3l in 1969 gave helicity its name and proved it is an invariant of ideal fluid flow. Arnol'd 
[HI proved several invariance results about helicity in 1974 (an English translation appeared 
in 1986). All three needed to invert the curl operator, or equivalently to compute d^a 
for a 2-form a, and utilized a version of the Biot-Savart operator. For more background 
on helicity and fluid dynamics, see for instance the book [|2l Chapter III] and the recent 
overview [fTTTl . 
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Helicity for vector fields is analogous to the linking number of a link and to the writhing 
number of a knot; many authors have explored these connections; see for instance 01 |2j 
[3l [5l [T2l [T4l|. For a cohomological view of helicity, involving the Biot-Savart operator for 
(k + l)-forms on (2k + 1) -dimensional subdomains of Euclidean space, see (SI. 

In |[T0l . DeTurck and Gluck extend the Biot-Savart operator to the three-dimensional 
sphere S 3 and hyperbolic space H 3 . We continue the investigation of how this story 
changes in the presence of curvature by looking at perhaps the most natural analog of 
the Euclidean setting, subdomains O C S 3 . 



1.1. Our aims. In this work, we develop an approach to electrodynamics on bounded 
subdomains of the three-sphere via the Biot-Savart operator. We provide integral formulas 
for Maxwell's equations and show precisely when BS acts as a right inverse operator to 
curl. Our formulas are geometrically meaningful, in that their integrands are preserved by 
orientation-preserving isometries of S 3 . Though verifying that Maxwell's equations hold 
on orientable 3-manifolds is an elementary exercise in differential forms, there exists a 
paucity of geometric formulas for electrodynamics in the presence of curvature; we aim to 
remedy this situation. 

Electrodynamics on the entire 3-sphere was developed in ifTOl . with applications to geo- 
metric knot theory and to the helicity of vector fields. Electrodynamics on compact subdo- 
mains (with boundary) of S 3 is the natural analog of the Euclidean setting, and it raises a 
rich and interesting set of issues: 

• The Hodge Decomposition Theorem for vector fields on is more complicated 
than on the three-sphere, because curl is no longer a self-adjoint operator and di- 
vergence is no longer the (negative) adjoint of gradient. 

• Current flows on bounded domains can deposit electric charge on boundaries and 
thereby affect Maxwell's equations. 

• Nonsingular current flows can be restricted to tubular neighborhoods of knots, en- 
abling connections between the writhing number of the core knot and both the 
helicity and flux of these flows. 

For modeling helicity and MHD on open, unbounded domains, such as the exterior of 
the sun in a Euclidean universe, the methods referenced above have had limited success. 
We propose an alternative for solar physics: modeling the universe as a giant three-sphere, 
and viewing the exterior of the sun as S 3 minus a ball, which forms a compact domain. The 
results herein might produce a useful model of solar plasmas. 



1.2. Listing of main results. The next section surveys the known results for the Biot- 
Savart operator in M 3 and in S 3 . Then, we demonstrate in Theorem 13.11 that Maxwell's 
equations hold on subdomains of S 3 via our setting, where we take a current V, its in- 
duced magnetic field BS(V) and the appropriate electric field. We then generalize one of 
Maxwell's equations: 
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Theorem 13.31 For Q a compact subset ofS 3 with smooth boundary, and V a smooth vector 
field on Q, then 



V, x BS(V)(y) 



V{y) inside f2 
outside £1 



— V y / 0o (Vj; ■ V(x)) dx + V y / 0o V"(ar) ■ n d(area x ). 
Jn Jan 



This theorem implies that for V divergence-free and tangent to the boundary, the curl 
operator acts as a left inverse to the Biot-Savart operator. These vector fields are the only 
ones for which BS inverts curl, as we show in Proposition 14.71 Any such V in this space 
that is also an eigenfield of BS must furthermore be an eigenfield of curl, with reciprocal 
eigenvalue, i.e., if BS(V) = XV, then Vx^ = (1/A)V. 

Also in section|4j we show that BS is a self-adjoint operator. We describe its image and 
find its kernel. 

Theorem 14.21 The kernel of the Biot-Savart operator on Q is precisely the subspace of 
gradients that are always orthogonal to the boundary dQ. 



2. Background on the Biot-Savart operator 

On October 30, 1820, Jean-Baptiste Biot (1774-1862) and his junior colleague Felix 
Savart (1791-1841) [|4l announced their landmark experiment which deduced the magnetic 
field associated to a current flow in a vertical wire and gave rise to the law (OQ), which bears 
their names. Soon after, Ampere deduced that for steady currents, VxB = p, J, a law 
named for him. (Henceforth, we shall choose units so that both /i and e each equal 1.) 
For an excellent history of the origins of electrodynamics, see |[T8l . 



2.1. Euclidean results from the Biot-Savart operator. The Biot-Savart law (OQ) holds for 
currents J, which we regard as divergence-free fields that are tangent to the boundary of 
the containment region; i.e., these are fields with no gradient component, by the Hodge 
Decomposition Theorem for vector fields (see Q). The Biot-Savart operator extends the 
law to act upon all vector fields; for x E fi, a compact subdomain of IR 3 , define it as 

(2) BS(V)(y) = -j- I V{x) x V - X dx = [ V(x) x V0(x,y) dx, 

^ n Jn \y ~ x \ Jn 

where the function <f>(x, y) = 1/ (47r|y — x\) is the fundamental solution to the Laplacian. 

This operator was introduced in O, which showed that BS is a bounded, compact, and 
self-adjoint operator. Furthermore, they showed that BS inverted curl for divergence-free 
vector fields tangent to the boundary. Our Theorems 13.31 and 14.21 are the direct analogs of 
their results, suitably adjusted in the presence of positive curvature. 
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2.2. Defining the Biot-Savart operator on S 3 . In this section, we define the Biot-Savart 
operator on the three-sphere. We follow the approach of DeTurck and Gluck IfTOl . 

Let VF(Cl) denote the space of smooth vector fields on fi. For a current J e VF(fl), its 
magnetic field is smooth except across the boundary d£l, where the field remains contin- 
uous; magnetic fields are linear in J. Therefore to define BS on the three-sphere we will 
work in the category of smooth, linear operators. 

To motivate the definition, let us first understand what essential properties a magnetic 
field in IR 3 possesses. First, magnetic fields are divergence-free. Second, Ampere's Law 
dictates that, for a steady current J, the curl of its associated magnetic field must equal J. 
Finally, we must describe how the Biot-Savart operator acts to gradients on S 3 ; for natural 
reasons, we require that it does so trivially. 

Let us explain this last requirement. In the Euclidean subdomain case, the kernel of 
BS is precisely the complement of the image of curl. On S 3 , the Hodge Decomposition 
Theorem for vector fields is straightforward; the image of curl is precisely the space of all 
gradients. Thus all gradients should lie in the kernel of BS. 

Theorem 2.1 (Hodge Decomposition Theorem on S 3 , cf. [19, Theorem 6.8]). Vector fields 
on S 3 decompose as 

VF{S 3 ) = K{S 3 )®G{S 3 ), 
and we call these spaces knots and gradients, respectively. Furthermore, 

K(S 3 ) = ker div = image curl 
G(S 3 ) = image grad = ker curl. 

We note that the space of steady currents is precisely the divergence-free fields, K(S 3 ). 
A special three-dimensional subspace is the set of left-invariant vector fields. These are all 
divergence-free, and are eigenfields of curl: V x U = — 2£7, for left-invariant U. 

The three properties above depict our definition. 

Definition 2.2. The Biot-Savart operator BS : VF(S 3 ) -> VF(S 3 ) is defined to be the 
smooth linear operator satisfying these three properties: 

1. BS is divergence free V ■ BS{V) = 

2. For V e K(S 3 ), curl inverts BS(V) V x BS(V) = V 

3. BS vanishes on gradients BS(Vf) = VV/ E VF(S 3 ) 

A straightforward exercise shows that BS is uniquely determined among smooth linear 
operators from VF(S 3 ) to VF(S 3 ). Suppose not, with B\ and B 2 two distinct operators. 
Then Bi(V) — B 2 (V) is both curl-free and divergence-free for all V. Hence by the Hodge 
Theorem, it is zero; thus B 1 = B 2 . 

Another approach defines BS using the Green's operator on vector fields, which inverts 
the Laplacian. The Green's operator commutes with all partial derivatives and hence with 
div, grad, and curl. See [fTBI for a description of Green's operators and Laplacians for vector 
fields of 3 -manifolds. 
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Proposition 2.3. The Biot-Savart operator is the negative of the curl of the Green's oper- 
ator: 

(3) BS(V) = -V x Gr(V) 

Proof. We show that this formula satisfies the three properties in Definition 12 .21 The oper- 
ator —V x Gr(V) is divergence-free because it is the curl of a vector field. Taking the curl 
of this operator yields 

V x (-V x Gr{V)) = A(Gr(V)) - V (GV(V • V)) 

(4) V x (-V x Gr(V)) = V - V (GV(V • V)) , 

where A is the Laplacian on vector fields (see f!6| for more details about this operator 
on three-manifolds). So if V is divergence-free, then the iterated curl above returns V, as 
desired, n.b., above we must also utilize the scalar Green's operator on functions. 

The last condition states that the operator should vanish on gradients. All gradients are 
mapped to zero because —V x Gr(Vf) = — GV(V x V/) = — Gr(0) = 0. □ 



2.3. Integral formulas for the Biot-Savart operator. The integral formula © for the 
Biot-Savart operator in Euclidean space requires the addition of vectors lying in different 
tangent spaces. To obtain an analogous formula on the three-sphere, we must decide how 
to move vectors among tangent spaces. Two natural choices exist: parallel transport along 
a minimal geodesic or left translation (or right translation) using the group structure of S 3 , 
which we may view as the group of unit quaternions or as SU(2). Each method has its 
advantages and disadvantages; we utilize both throughout this article. 

Let P yx V denote parallel transport of V E T X S 3 to the tangent space T y S 3 ; likewise let 
(L yx -i)*V or L*V denote left-translation. Let a(x, y) be the distance on the three-sphere 
between x and y. 

Regarding notation, it is important to keep track of at which tangent space a derivative is 
taken. We affix a lowercase subscript to the nabla operator to accomplish this, e.g., V x f is 
the gradient of / at x. 

Theorem 2.4 ( IITOl Theorem 2]). As an integral in which vectors are moved via parallel 
transport, the Biot-Savart operator is given by the formula: 

(5) BS(V)(y) = f P yx V(x)xV y( j ) (x,y)dx. 

Jn 

As an integral in which vectors are moved via left translation, the Biot-Savart operator 
is given by the formula: 



(6) 



BS(V)(y) = / (L yx -x)*V(x) x V„0 O dx 
In 



An 2 



(L yx -i)*V(x) dx + 2V y / (L yx -i)*V(x) ■ V y (j>i dx. 
n Jn 
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This theorem expresses BS(V) as a convolution of V, thought of as an electric current, 
with (the gradient of) appropriate potential functions. The left-translation version of BS 
requires two different convolutions, while the parallel transport version requires only one. 
The potential functions are 

<p(a(x, y)) = --^2 (7r - a) csc(a) 0oOO, y)) 

<Pi(a(x,y)) 

The function O is the fundamental solution of the Laplacian on S 3 and defines <pi as 
follows: 

A0 O = 5(a) - — A0i = 0o - [0o]- 

Here 5(a) represents the Dirac delta function and [/] denotes the average value of / on 
S 3 . Recall that the Laplacian of a function has average value zero on a closed manifold; 
hence we require the constant \/2n 2 . 

Likewise, the function <p is the fundamental solution of the shifted Laplacian (A — 1); 
i.e., A(f) — 4> = 5(a). We note for future reference that 4> (a) = 4>(a) cos a. 

Using the Biot-Savart operator, DeTurck and Gluck IflOl obtain explicit integrals for 
calculating the linking number of two knots on S 3 . They also generalize Calugareanu's 
formula Link = Twist + Writhe [|9l by finding explicit integrals for the twist and for 
the writhe of a knot on S 3 . 



3. Electrodynamics on subdomains of S* 3 

This section investigates electric and magnetic fields on subdomains of the three- 
sphere via the Biot-Savart operator. We assume f2 is compact with piecewise smooth 
boundary. We define BS to be the same operator from Definition |2.2[ with integral for- 
mulas given by © and ©; however, the region of integration is now restricted to f2. We 
may then view this operator as BS : VF(tt) — > VF(Q). For y Q, these formulas also 
define the behavior of BS(V) outside of the domain. 

3.1. Maxwell's equations. We begin by considering a current V(x) on Q and computing 
BS(V) as its corresponding magnetic field. If the vector field V either has a nonzero 
divergence or is not tangent to the boundary <9f2, then it no longer represents a steady 
current contained in £1. In this case we take a time-varying electric field, so that our system 
is 'closed'. If V is not divergence-free, charge will accumulate at a rate of minus the 
divergence, so let p(x,t) = — (V ■ V)t represent the volume charge density in f2. If V 
flows across the boundary, charge will accumulate on the boundary; let a(x, t) = (V ■ h)t 
represent the surface charge density on dVL. Each of these charge distributions determines 



1 

471-2 
1 
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an electric field: 

(7) E p {y,t) = -\VvJ <!>o {V x ■ V{x)) dx) t = V y j fopdx 

(8) E a (y,t) = [V y / 0o V(x) ■ n d(area x ) J t = V y / (f) a d(area x ) 

We can view dE p / dt as the time rate of change of the electrodynamic field due to dp/ dt, 
the change in volume charge density; also we can view dE a /dt as the time rate of change 
of the electrodynamic field due to da/dt, the change in surface charge density. We also 
may view dE p / dt as an electrostatic field itself due to the time-independent volume charge 
dp/dt — V • V; similarly, we may view dE a /dt as an electrostatic field due to the time- 
independent surface charge dp/dt = V • ft. In this section we adopt the former viewpoint, 
while later we will make use of the latter. 

We consider the sum of these as the electric field E(y,t) = E p + E a . The following 
theorem demonstrates that this electrodynamic setup obeys Maxwell's Equations. 

Theorem 3.1. Let V = V(x) be a smooth vector field on Q; we view V as a steady 
current. For the electric field E(y, t) defined above and the magnetic field B{y) = BS(V) 
corresponding to V, Maxwell's four equations hold at each y G Q: 

1. V-E = p 3. V-5 = 

dB dE 

2. V x E = — 4. VxB = y+- 

dt dt 

(n.b., we have chosen units so that both the permittivity of free space e and the permeability 
of free space /i are identically 1.) 

By the end of this subsection, we will have proven Theorem l3.ll The following proposi- 
tion restates Maxwell's third equation. 

Proposition 3.2. For V a smooth vector field, BS(V) is divergence-free, whether defined 
on Q or on S 3 . 

Though Definition [2T2] demands this to be true for BS, we state it here since we provide 
an independent proof of this fact for its parallel transport formula © in section [3~4l 

Next, we compute the curl of BS: 

Theorem 3.3. For Q a compact subset ofS 3 with smooth boundary, and V a smooth vector 
field on Q, then 

v t x Bsmw = { v n {y) im i de f 

^ (J outside i I 



-V y 0o (V-,. ■ V(x)) dx + Vy 

Ju Jd 



0o V(x) ■ h d(area x ). 

an 
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If V is divergence-free and tangent to dQ, notice that curl acts as a left inverse to BS, 
(10) V, x BS(V)(y) 



V(y) inside 
outside f2 



The above theorem is useful in solving certain energy minimization problems for vector 
fields on S 3 (e.g., analogs of the Woltjer and Taylor problems in plasma physics). Solutions 
are often eigenfields of curl and, by (flOl) . also of BS (assuming they have no gradient 
component). 

We provide two proofs of this theorem in section 13 .61 one using parallel transport of 
vectors and one using left-translation. The first argument relies on a result, the Key Lemma, 
described in section [331 

Recognizing the last two terms of © as dE p /dt and dE a /dt, respectively, we obtain as 
an immediate corollary that Maxwell's fourth equation holds. 

Corollary 3.4. For y £ Q, 

dE dE 

(11) VyXB S(V)(y) = V(y) + —t+ " 



dt dt 

Remark 3.5. When f2 = S 3 , Theorem 13.31 above reduces to 

(12) V, x BS(V)(y) = V(y) -V y I O (V, • V{x)) dx. 

In this case Proposition |2.3| and (HJ) provide an easy proof: 

(13) VxBS(V) = V x (-V x Gr(V)) = V -VGr(V -V). 



3 



The last term is precisely the gradient in (PT21) . Alas, the Green's operator for Q ^ S 
can be much more complicated, and we will require a more difficult argument to prove 
Theorem [33] in general. O 

Example 3.6. There exists a three-dimensional subspace of VF(S 3 ) consisting of all left- 
invariant vector fields. Any such field U is divergence free and satisfies V x U = —2U. On 
S 3 , we calculate that BS(U) = — (1/2)Z7, and so V x BS(U) = U, as expected by CU). 

Consider the solid torus Q a that consists of all points within spherical distance a of the 
circle x 2 + y 2 = 1, which we denote as 7. Let u\ be the unit, left-invariant field with an 
orbit along 7. Then, BS(u\) = — \u\ + \ cos 2 a W, where W is the unique curl-free, 
divergence-free field that has unit length along 7; all orbits of W are longitudes of £l a . (In 
terms of the Hodge Decomposition Theorem (section I47TT) , W is a harmonic knot.) Then, 

j cos 2 a 

V x BS (ui) = --Vxm! + — ^— V x W = ui, 

as expected by Theorem 13 . 3 1 since ui is divergence-free and tangent to dfl a . O 



We now return to prove Theorem l3.ll 
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Proof of Theorem \3J\ We have already seen through Proposition |3.2| and Corollary [3]4] that 
the last two of Maxwell's equations hold in our setting. We must now show the first and 
second equations do as well. 

The divergence of E in Vt can be calculated by taking the divergence of both sides of 
Maxwell's fourth equation: 

V-VxBS(V) = V-V + V- — 

at 

o = v-v + ^iy-E) 

Since the current V does not depend upon time, we integrate to obtain that V • E = — (V • 
V)t = p, Maxwell's first equation. 

Maxwell's second equation holds trivially. The electric field E is defined by two gra- 
dients and ergo its curl is zero. The magnetic field BS does not depend upon time, even 
when the electric field E is time-dependent; hence V • E — = dB/dt. □ 

The remainder of section [3] is devoted to proving the results claimed above. First, we 
need to define a computational aid in R A , the triple product. 

3.2. Triple products in R 4 . Let A, B, C be vectors (or vector fields) on R 4 . Let a e [0, 7r] 
be the angle between vectors A and B. 



Definition 3.7. The triple product of A, B, and C is the vector in 

[A, B,C}= det 



d\ 0,2 CL3 0-4 

h b 2 63 h 

Cl c 2 c 3 c 4 

Xi X 2 £3 ^4 



The triple product of three vectors in R 4 is the analog of the cross product of two vectors 
in R 3 . Indeed, the product of n — 1 vectors in R n is similarly defined as the determinant of 
antixn matrix. 

Three useful properties of triple products are that 

(1) They are multilinear and alternating, i.e., [A, B, C] = [B, C, A] = -[A, C, B}. 

(2) [A, B, C] is orthogonal to A, B, and C. If A, B, and C are linearly independent, 
then {A, B, C, [A, B, C] } forms a basis that agrees with the standard orientation on 
R 4 . If A, B, and C are linearly dependent, then [A, B, C] = 0. 

(3) If A is a point in S 3 (i.e., \A\ = 1), and B and C are tangent to the three-sphere at 
A, i.e., B,C G TaS 3 , then [A, B, C] lies in TaS 3 and is equal to the cross product 
B x C. More generally, for A e S 3 , the triple product [A, B, C] = B L x C L , 
where B 1 - = B — (A ■ B) A E TaS 3 , is the component of B perpendicular to A 
(and likewise C 1 - e T A S 3 ). 

Often, calculations require a formula for an iterated double product. 
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Lemma 3.8. Let A, B, C be vectors in R 4 . Let C 1 - represent the component of C that is 
orthogonal to the plane spanned by A and B. 

(14) [A,B, [A,B,C]] = -\A\ 2 \B\ 2 sm 2 (a) C L 

Proof. Let a be the angle between A and B. We use Gram-Schmidt to find B 1 - orthogonal 
to A and to find C 1 - orthogonal to both A and B^: 

(15) B ± = Bsin(a) 

\B\(A ■ C)-\A\ cos a(B-C) \A\{B ■ C) - \B\ cos a{A ■ C) 

(16) G = G |, mnl , ■ 2/ \ A iTiTTTTo • 2/ x B 



Let D = [A, B, C]. Assume {A, B, C} are linearly independent, else D — 0. Then D 
is orthogonal to the span of A, B, C and the basis {A, B, C, D} has positive orientation in 
R 4 . The length of D is \D\ = \A\ \B\ \C^\ sina. 

Let E = [A, B, D] = [A, B, [A, B, C]]. Then E is orthogonal to D, so it is a linear 
combination of A, B, and C. Since E is also orthogonal to A and B, it must be a multiple 
of C L . To ensure that the basis {A, B, D, E} has positive orientation in IR 4 , the vector 
E must point in the direction of — C L . The length of E is \E\ = \A\ \B\ \D\ sin a = 
|A| 2 |_B| 2 \C ± \ sin 2 a, which proves the lemma. □ 

3.3. The calculus of parallel transport. Let x, y be non- antipodal points on S 3 ; we will 
view them as unit vectors in R 4 . Recall, a is the distance on S* 3 between them. The 
component of y perpendicular to x is y 1 - = (y — x cos a). 

Since the three-sphere has unit radius, the gradient of a with respect to x (in T X S 3 ) must 
be a unit vector which points away from y along — y 1 . Thus 

nn\ V7 I \ xcosa-y 

(17) v x a{x,y) = 



(18) V y a(x,y) 



1 2/ | sma 
—x 1 - y cos a — x 



\x^\ sma 



For any two unit vectors x,y £ R n such that x ^ ±y, the unique map M E SO(n) that 
maps x to y and fixes all vectors orthogonal to both x and y is 

M(V)=V ( V -( x + y)) xl (V-x)(l + 2(x-y))-(V-y) 
V ' l + (x-y) l + (x-y) 

For V a tangent vector at x E S 3 , this map M precisely describes its parallel transport 
to the tangent space at y. The expression above simplifies to 

(19) P yx V = V- i { ^ {x y) y) (x + y) = M(V) 

Remark 3.9. For a vector vatxES 3 that points parallel to the geodesic 7 running through 
x and y E S 3 , left-translation from x to y is exactly the same as parallel transport. The two 
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methods differ only in how they treat components that are perpendicular to the geodesic 7. 
Hence we compute 

(20) V y a = -P yx V x a = -(L yx -i)*V x a. 

The interested reader is invited to show this directly using the group structure of S 3 and 
equations CH]) and (fT8T> . O 

Remark 3.10. For any function f(a) that depends only upon the distance a(x, y) from x 
to y, its gradient with respect to x variables is V x f(a) = f'(a)V x a. Thus the methods of 
transporting its gradient vector are equivalent, 

V y f(a) = -P yx V x f(a) = -(L yx ^V x f(a). O 



3.4. The divergence of BS. Though Definition 12.21 requires the Biot-Savart operator to 
have zero divergence, here we furnish an independent proof of this fact; we utilize the 
parallel transport formula © for BS. 

Proof of Proposition \3. 21 We start by taking the divergence of ©, 

V y -BS{V){y) = Vy f P yx V x V y <f> dx = [ % • {P yx V x V v <j>) dx 

Jn Jn 

Apply the vector identity V ■ (A x B) = (V x A) ■ B - ( V x B) ■ A to the integrand: 

(21) V y ■ BS(V)(y) = fv v <p- (V„ x P yx V) dx, 

Jn 



since the curl of V y 4> is zero. Next we compute the curl of P yx V using (|T9l) : 

V y x P yx V = V y x V(x) -V y x - V ' V (x + y) 

1 + cos a 

Since x is fixed with respect to y, V y x x = 0. Applying the Leibniz rule for curls, we 
obtain 

/ V ■ y \ V ■ v 

V y x P yx V = -V y x (x + y) - V y x (x + y) 

\ 1 + cos a J 1 + cos a 

The cross product above is taken in the tangent space at y, so a vector crossed with y 
contributes nothing. Similarly V y x y = 0. 

(22) V 9 x P SI V {X ) = -V, (j^) x x 

Now calculate the gradient above: 



y ^l + cosa/ 1 + cosa (1 + cosa^ 2 y 



V y (y • y) + -— — V y cos a 
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To proceed, we need to calculate Vy (V ■ y). It should point along the component of 
V that lies perpendicular to y. Let 9 be the angle between V(x) and y in R 4 . Since 
( V ■ y) = \V\ cos 9, we calculate 

Vy(y-y) = ^(|r|cos0)^ = -|V|sm^ COS ^ 1/|M 



sin 9 



(23) V y (V-y) = V-{V-y)y = V~ 



A similar formula holds for the gradient of the inner product of any two vectors in R™. As 
a check, we compute cos a = — smaV y a = x L . by (fT8l). which agrees with (1231). 
Then, 



"\l + cosa/ 1 + cosa (l + cosa) J 



After substituting this in, we may express equation (|22l) as a triple product, described in 
section l3~!2l 

V, x P yx V = ( -— ^ ^ + ^ Nn x x 



Vj X Py X V 



1 + cos a (1 + cos a) 2 



x xx 



1 + cos a (1 + cos a) 



A triple product is zero if the three vectors are not linearly independent, so the x L term 
above vanishes. 

V y x P yx V = — ^ [y,V,x] 

1 + cos a 

Recall that the triple product is unchanged if we switch to vectors perpendicular to y: 
[y, V, x] = [y, V ± , x 1 -]. Again, let the triple product represent a cross product. 

V, x P yx V = -— L_ [y, V\ x 1 ] = E T y S* 

1 + cos a L J 1 + cos a 

We finally return to (T2~Ti 

V, • BS(V)(y) = ( V y <t>- (Vy x P yx ) dx= f <j)'(a) V y a ■ ( V± t * X± ) dx 
Jn ' Jn \l + cosaJ 

(V 1 x x 1 -) dx, 



Via) —x~ 



1 + cos a \x J 

where we have substituted (fT8l) . The product x 1 - ■ (V 1 - x ar 1 ) = 0, so the integrand vanishes 
identically at each point. We conclude that the Biot-Savart operator is divergence-free. □ 
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3.5. Key Lemma. In this section, we prove an important lemma relating vector fields and 
functions on S 3 . DeTurck and Gluck IfTOl have also proven this lemma; we furnish an 
independent proof. 

Lemma 3.11. [Key Lemma] Let x,y be two non-antipodal points in S 3 C M 4 . Let 4> = 
4>(a) be a function, depending only on a, which may have a singularity at a = but is 
otherwise smooth. Let V be a tangent vector at x G S* 3 . Then, 

(24) V y x {P yx V x V,0} - V y {V ■ V x (0cosa)} = (A0 — (f>) (V — (V ■ y) y) . 

Here the Laplacian is taken in terms functions on S 3 , where a can be taken as a coor- 
dinate. We will use the Key Lemma in the next section to prove Ampere's Law. Another 
application of the Key Lemma is in proving an integral formula for the linking number of 
two knots on S 3 ; refer to JlOJ for details. 



Proof. Denote the two terms on the left-hand side of the Key Lemma as terms C and G, 
respectively. We will show that C — G = (A<fi — <p) (V — (V ■ y) y). 

We convert the cross product in C into a triple product and utilize Remark [3 .101 to com- 
pute the gradient of 0, 



(25) 



C = Vy X 



y, P yx V, -r— (ycosa - x) 



sm a 



We calculate the parallel transport of V(x) via (fT9l) , and the triple product above becomes 

V-y 



[y, PyxV, (ycosa- x)] 



y, v 



1 + cos a 



[x + y) 



-x 



[y, v, -x] 



since the triple product is linear and the terms with repeated variables vanish. Hence, C 
becomes 

C = V„X ^[y, X , V ], 



sm a 



Applying the Leibniz rule for curls, we obtain 

>'(«) 



(26) 



C = V, 



sm a 



x [y, x, V] + V y x [y, x, V] 



sm a 



Call these two terms C\ and C 2 respectively, so C = C x + C 2 . We analyze them sepa- 
rately. The first one, C\, is converted to an iterated triple product. 



v^v^Ay^.v] 



sin a 



d 



da V sm a 



y cos a — x r . 
3/> > x ' V \ 



C x 



L >" sin a — 4>' cos a 
sin 3 a 



sm a 
[y, x, [y, x, V}} 
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Now we utilize Lemma [3781 to evaluate the iterated triple product. In order to do so, we 
must first calculate V*, the component of V orthogonal to x and y from equation ( TT61 ). 



Then G\ becomes 



V = V - y + — — (V ■ y) x 

sin a sm a 



k " sin a — d>' cos a 



Ci = 5 -sin a V ) 

sin a 



Jl sin a — d>' cos a 



(27) d = y 3 ^ (sm 2 aV+(V-y)y-cosa(V-y)x) 



sin a 



Now we consider C 2 from (|26l) . Since [y, x, V] = [x, V, y] = jx- 1 , V*, y], we can rewrite 
it as 

(28) C 2 = ^V y x[y,x,V] = ^V y x[x\V\y}. 

ana sina 



Claim. V y x [x ± , V*, y] = 2(x ■ y)V - 2(V ■ y)y. 

The claim follows in a bare-hands computation similar to the ones above. Start by ex- 
pressing the triple product as a curl, and then use vector identities and (|23l ) to compute all 
necessary derivatives. We leave this to the reader. 

Applying the claim to (|28l) . we obtain 

C 2 = 2cj)' C —V-2(j)'^—{V-y)x 
sm a sm a 



Summing terms C\ and C%, we obtain the following expression for C: 

c = U" + ^) v 

V sin a / 

\ sm a sm a sin a / 



Next, we turn our attention to the second term G — V y (V ■ V x (<fi cos a) ) in the statement 
of the Key Lemma. 

( d \ ( x cos ot — y \ 

G = — V y ( V ■ — ((pcosa) V x a ) = V y ( (0'cosa — 0sina) V ■ 



da J V sin a 
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Since V E T X S 3 , the product V - x = 0. Thus, 



G = ~Vy 



cos a — fflsma , Tr . 

: (V -y) 

sin a 



„ , TJ . N „ / ..cos a ,\ /..cos a , \ __ /T s 

G = -(V-y)V y U'- 0-0- <P)V y (V-y) 

V sin a / V sma / 

„ , , d ( , cos a \ / , cos a \ 

(30) G = -(V-y)-U'- Vj/ a- 0'- 0)V tf (y.y 

aa V sma / V sma / 



By (|23l . V w (y • y) = V - (V • y) y. We substitute via (QH) for V y a to obtain 
(3D 



G = - (V • y) ( 0"— ^- - 0'^- + 0' ) (ycosa - x) 

sin a sin a sin a 



: (<//—-</>) (^-(V.y)y) 
sm a / 



Finally, C — G, which forms the left-hand side of the Key Lemma statement, is merely 

the difference of equations (|29l) and (|3TT) . All x terms cancel, leaving 

,cosa \ / ,// ,,cosa \ , 

'AV + -0"-20'- + ){V-y y 

/ V sma / 



C-G 



sm a 



(32) C-G = U" + 20'^ - cf>) (V - (V ■ y) y) 

V sm a / 

cos a 

In terms of the coordinate a, the Laplacian of 0(a) on S 3 is A0(a) = 0" + 26' . Thus 

sma 

we have proven the Key Lemma, since 

C-G = (A0-0)(F-(V-y)y) . □ 

3.6. Proving Theorem 1331 

3.6.1. Parallel transport proof of Theorem [Ol The Key Lemma, as perhaps its most im- 
portant consequence, directly proves Theorem |3.3l Maxwell's fourth equation, which states 



V, x BS(V)(y) 



V(y) inside 
outside 



in Jan 



y(x) • n d(area x ) 



Proof. We begin by inserting the function 0(a) = — ^(tt — a) esc a, which appears in the 
parallel transport formula © for BS, into the Key Lemma. Recall that 0o(a) = 0(a) cos a 
and A0 — = 5(a) = 6(x, y). The Key Lemma states 

m V y x {P^y x V,0} - V„{V ■ V, (0cosa)} = (A0 - 0) (V - (V ■ y) y) . 
Now, integrate both sides over f2 with respect to x: 



[ V y x {P^V x V y 0} tfe - / % {V ■ V X O } dar = / , 



S(x,y) {V{x) - (V -y) y) dx 
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We may interchange the integral in x variables with the gradient and curl operators on 
the left-hand side, since they are in terms of y. 

V y x [ P yx V{x) x V y (j) dx-V y [ V (x) ■ V^ dx = [ 5(x, y) (V(x) — (V ■ y) y) dx 
Jn Jn Jn 

The first term on the left-hand side is precisely the curl of the Biot-Savart operator. Let 
us substitute that into the equation above. 



(33) V,, x BS(V) =V y [ V ■ V x (p dx + [ 

Jn Jn 



y) (V(x) - (V(x) ■ y) y) dx 



We now apply the vector identity V ■ V0o = V • ((f>oV) — </> V ■ V to expand the first 
RHS integral: 

/ V ■ V x <p dx = I V x ■ (0 o y) dx - Vy I 0o V x • V dx 
Jn Jn Jn 

4>o V ■ n d(area x ) — V y / 0o V x • V dx 



'an Jn 
after applying the Divergence Theorem. Substituting into (|33l , we have found three of our 
four desired terms. 

V y x BS(V) = V y 4> V ■ hd(area x ) - V y / (p (V 2 ■ V) dx 
Jon Jn 

+ / 8(x, y) {V(x) - (V(x) ■ y) y) dx 



Jn 

To evaluate the last integral, we substitute y for x in its integrand to obtain V(y) — (V(y) ■ 
y) y. Since V(y) lies in the tangent space of y, this reduces to just V(y), which is assumed 
to be zero outside of £1. We write this fact explicitly, 

S( X , y) (V( X ) - . y) y) dx = { V ^ , 

which completes the proof of the theorem. □ 



3.6.2. Left-translation proof of Theorem LOl 

Proof. The left-translation formula for BS is comprised of three integrals: 

BS{V){y) = f {L yx -x)*V{x) x V y o dx 
Jn 

^ J {L yx -i)*V(x) dx + 2V J/ J (L yx -i)*V(x) ■ V^i dx 



(6) 



We take the curl of each term separately; the first integral requires significant analysis 
while the others do not. The third is a gradient; its curl must be zero. Meanwhile, the 
second term produces a vector field — \\V\, where [V] measures the average value of V 
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on all of S 3 ; such a vector field is left-invariant. Recall that all left-invariant fields are 
eigenfields of curl, so V x — \\V\ = [V]. 

Now we calculate the curl of the first integral, which for convenience we denote as B. 
We truncate notation to L*V in lieu of (L yx -i)*V(x). 



V y x B = / Vy x (L*V x V y (f) ) dx 
Jn 



Apply the vector identity for the curl of a cross product: 

Vx(U xW) = [W, U] + (V • W)U - (V ■ U)W, 
where we set U = L*V and W = V y 4>Q. We obtain three terms to analyze. 

(34) V y x B = [ [W,U]dx + [ (A y o ) L*V dx - f (V„ ■ L*V) V y o dx 
Jn Jn Jn 



The third term of (1341) . We claim this term is zero. For fixed x and varying y, the vector 
field L(y) = (L yx -i)*V(x) is a left-invariant vector field. Left-invariant vector fields on S 3 
are divergence-free, so this integrand of this term vanishes pointwise. 



The second term of (1341) . The function O was defined to be the fundamental solution of 
the Laplacian, i.e., A(p (a) = 5(a) — l/2n 2 . Then, the second term can be rewritten as 



I (A y (f) ) L*V dvol x = I 5(a) (L yx -i)*V(x) dx - — f (L yx 



-i)^V(x) dx 



Since 8(a) vanishes except when a = 0, i.e., when x = y, the first RHS integral simply 
takes on the value of the vector field when x = y. The second RHS integral is the formula 
for [V] . Thus the second term of (l34l) is 

(35, / n ( W ,V^(^> ™ )- lV] 

The first term of (|34|) . We will need this lemma, proven in ifTTl . 

Lemma 3.12. On S 3 , let U be a left-invariant vector field, let W be any smooth vector 
field, and let G be a gradient. Thenu 

V W U = W xU and V(U ■ G) = [U, G}. 

The lemma allows us to rewrite the first term as 

(36) f [W, U)dx = - I V y (W -U)dx = V y [ V^ • UV dx 

Jn Jn Jn 



Regarding notation, we utilize uppercase subscripts on V to denote covariant derivatives, e.g., Vy/U . 
Lowercase subscripts, e.g. V y , continue to refer to the point at which the nabla operator is applied. 
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Recall from (|20l) that the gradient of <fr changes sign when switching variables from y 
to x appropriately, V^0o — ~ (Ly X -i)*Vx<j>o- Therefore, 

/ [W,U]dx = -V y I -{L yx -^V x (t> G -{L yx -i)*V{x)dx 
Jn Jn 

= V y V x <po ■ V{x) dx, 
Jn 

since there is no need to left translate before taking the dot product. Now, we apply the 
Leibniz rule for divergences to obtain 

/ \W, U] dx = V y [ V x ■ <j) V(x) - O (V, • V(x)) dx 
Jn Jn 

(37) / [W,U]dx = V y <f) iy x -V) dx - V y i <p V ■ n d{area x ). 
Jn Jn Jan 

Gathering our results from (1351) and (1371) . along with the [V] contribution from the second 
integral that comprises BS, we have proved Theorem [331 

^ x * = ( T S ) - M 

+ V y / 0o (Vc -V) dx — Vy 4>qV ■ n d(area x ) 
Jn Jan 

%x-I[V] = [V] 

The sum of these two terms is the curl of BS, as desired. □ 



4. Properties of Biot-Savart on subdomains 



In this section we calculate the kernel of the Biot-Savart operator and prove a result 
about its image. As an aid to applications, we determine for which vector fields is curl a 
left-inverse to BS. Finally, BS is self-adjoint, bounded operator. We prove self-adjointness 
in section |44l our upcoming paper |fl"5l proves boundedness. 

To discuss the kernel and image, we first need to understand how vector fields decompose 
on subdomains of the three-sphere. 

4.1. Hodge Decomposition Theorem for c S* 3 . In order to discuss the properties of 
BS, we require an understanding of vector fields on subdomains of S 3 decompose. The 
Hodge Decomposition Theorem provides this, in relation to the operators div, grad, and 
curl. See for a detailed proof for vector fields in R 3 ; slight modifications to that argu- 
ment prove the following version on S 3 . 

Theorem 4.1 (Hodge Decomposition Theorem for S 3 ). Let be a compact, three dimen- 
sional submanifold of S 3 with dVL piecewise smooth. Then, there exists a decomposition of 
VF(fl) into five mutually orthogonal subspaces, 

VF(Q) = FK © HK © CG © HG © GG, 

where, 
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FK 


= fluxless knots 


= {v 


■ V = 


0, V-n 


= 0, all interior fluxes = 0} 


HK 


= harmonic knots 


= {v 


■ V = 


0, V -n 


= 0, V x V = 0} 


CG 


= curly gradients 


= {v 


= V<p 


,v-v = 


= 0, all boundary fluxes = 0} 


HG 


= harmonic gradients 


= {v 


= V0 


,V-V-- 


= 0, <p locally constant on dfl} 


GG 


= grounded gradients 


= {v 


= V0 


, (Plan = 


0} 



The subspaces HK and HG are finite dimensional and 

HK Hl(tt,R) = R^nnsdn 

JJQ ~ JJ ~ jg^ | components of 9fi|— | components of f2| 

Furthermore, 

kerdiv = FK © HK © CG © HG 

image curl = FK © HK © CG 

kercurl = HK © CG © HG © GG 

image grad = CG © HG © GG 

We refer to the subspace K = FK®HK as the fluid knots, or simply as knots; it consists 
of all vector fields that are both divergence-free and tangent to the boundary. 

4.2. Kernel of the Biot-Savart operator. By definition, the Biot-Savart operator on S* 3 
maps the subspace of gradients to zero. No knot on S 3 lies in the kernel of BS or else Am- 
pere's Law (Theorem [3]3]) would fail. Hence the kernel of Biot-Savart on the three-sphere 
is precisely the space of gradients. Gradients on S 3 all behave like grounded gradients 
found on a compact subset C S 3 . There the Hodge Decomposition Theorem for vector 
fields is more complicated, so we ponder, how do the other subspaces behave? What is the 
kernel of BS on ft? 

Theorem 4.2. Let Q C S 3 be as described above. The kernel of the Biot-Savart operator 
on Q is precisely those gradients that are orthogonal to the boundary, i.e., 

ker BS = HG(Q) ®GG(Q). 

When discussing the kernel, we must note carefully that as an operator BS maps into 
VF(£l). Though we often extend BS to defining a vector field on S 3 — f2, that is not its 
natural target space. Therefore, a vector field V lies in the kernel if and only if BS(V) = 
inside f2; a priori nothing is known of its behavior on S 3 — f2. In proving Theorem 14 .21 
we will show that if BS(V) = throughout fi, then BS'tV) must vanish identically on the 
entire three- sphere. 

In order to prove this theorem, we require a few preliminary results. 

Lemma 4.3. Let Q be as above, and let n be the outward pointing normal vector on dQ. 
Consider a vector field V G VF(Q) and let y G S 3 . Then, 

/ (L y:r -i)* (V x x V(x)) + 2(L yx -i)*V(x) dx = - (L yx ~i)* (V(x) x h) d(area x ) 
Jn Jau 
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Proof. Let U(x) be any left-invariant vector field on S 3 . Apply the divergence theorem to 

V x U: 

(38) / V • (V xU) dx = I (VxU)-n d(area x ) 

Jn Jan 

Now examine the right-hand side of this equation. 

(V x U) ■ n d(area x ) = — I U(x) ■ (V(x) x h) d(area x ) 
uv. Jan 

{L yx -i)JJ{x) ■ (L yx -i)* (V(x) x n) d(area x ) 

an 

U(y) ■ (L^-i)* (V(x) x h) d(area x ) 

an 

-U(y) ■ / (L yx -i)* (V(x) x n) d(area x ) 
Jan 



Now examine the left-hand side of equation (1381) . 

I V • (V X U) dx = [ U{x) ■ V x V(x) - V(x) ■ V x U(x) d{area x ) 
Jn Jan 



Since U is left-invariant, V x U — —2U. 

dx 



I V • {V x U) dx = I U ■ (V x V + 2V) 
Jn Jn 

{L yx -i)JJ{x) ■ {L yx -i)* (V x V{x) + 2V{x)) dx 
U(y) ■ (V x V(x) + 2V(x)) dx 



U(y) ■ / (V-i)* (V x V(x) + 2V{x)) dx 
In 



The two sides of equation (1381) are equal, which implies 



U{y)- (I r i) t (VxV + 2V) dx = -U(y)- (L yx -i)*(V(x) x h) d(area x ). 
Jn Jan 

Since this holds for any left-invariant field U, we may conclude that the projections of the 
two integrals onto the space of left-invariant vector fields must be equal. Both are integrals 
of left-translated fields and hence both integrals are left-invariant vector fields depending 
upon y. Therefore, these two integrals must be equal, which proves the lemma. □ 

We require another lemma which involves a particular energy estimate. In section [3l we 
defined the electrostatic field E a due to the time-independent surface charge a = V ■ h 
on dVt. (We also viewed E a as the time derivative of an electrodynamic field © due to a 
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surface charge (V ■ h) t but do not adopt that view in this section.) For convenience, we 
now drop the derivative notation and hereafter refer to this field as 



(39) E a (y) =V y 0o V(x) ■ h d(area x ). 

Jan 

Call ip the potential function for the electrostatic field, E a = —Vip. 

Lemma 4.4. For V = V/ G VF(Q), let E a the electrostatic field that it generates as 
described above. Then the energy of E a is related to its potential ip as 

/ \E a \ 2 dy = j ip{y)Vf ■ n d(area y ). 
Js 3 Jan 

Proof. Let V = V/, and let a = V/ ■ n, which we regard as a surface charge on <9f2. Then 
we can write 



E a = V y / 0o cr d(area 2 
Jan 



Now, smoothly extend the surface charge o on dfl to a volume charge p on a thickened, 
compact neighborhood N{dVt) of the boundary dVt. We choose p to be C°° smooth with 
support in N(dCl). The electrostatic field E p resulting from p approximates E a and is 
expressed as 



(40) E p = V y (P pdx. 

JN(an) 

If we regard E p = —Vip p , then the divergence of E p isV-E p = — Aip p . We may extend 
the domain of integration of ip p to be all of S 3 , since p vanishes outside of N(dCl). 



Recall that (fi is the fundamental solution to the Laplacian on S 3 , hence Aip p = — (p — 
Thus, the divergence V • E p — p — [p] . 

We consider now the integral of ip times p, and again extend the domain of integration 
tobeS 3 : 



/ ip p dx = / ip p dx = ipV ■ E p dx + / ip [p] dx 
JN(an) Js s ' Js s ' Js a 

•dx = / V ■ (ip E p ) - Vip ■ E p dx + [p] tp dx 
Js 3 Js 3 



'N(an) 

(41) / ippdx = 0+ / \E p \ 2 dx + [p][<0] 

JN(an) Js 3 

Now we shrink the neighborhood N(dCl) so that it approaches the boundary d£l. Then 
p — ► a and [p] — > 0. Thus the LHS of (flTI) converges to f gn ip a d(area x ) while the RHS 
converges to f„ 3 \E a \ 2 dx. Thus we obtain our desired result, as equation (14TT) converges to 

ip a d(area x ) = / \E a \ 2 dx. □ 
an Js 3 
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One more result, the following energy estimate, is required before beginning the proof 
of the kernel of BS. 

Proposition 4.5. Let V be a divergence-free vector field on Q C S 3 , and let E a be its 

associated electrostatic field. Then, 



\E a \ 2 dy< / \V\ 2 dy 



Proof. When V is divergence-free and tangent to the boundary, the electrostatic field E a 
0. So it suffices to prove the proposition for V a divergence-free gradient, i.e., V E CG 
HG. We write V = V f , which implies / is harmonic. The preceding lemma states 



\E a \ 2 dy 



S3 



ip(y) V/ ■ n d(area y 



mi 



Apply Green's first identity: 



(42) 



\E a \ 2 dy 



S 3 



W ■ V/ + ipAf dy 



-E a -Vfdy 



This is the L? inner product of —E a and V/; apply the Cauchy-Schwarz inequality: 



-E a -Vfdy < 



\E„ 



< 



\E, 



?dy 
r \ 2 dy 



1/2 



iV2 r 



|V/| 2 dy 



S3 



\Vf\ 2 dy 



1/2 



nl/2 



Substitute this inequality into equation (|42)) . to conclude 



\Ea\ 2 dy < 



|K| 2 dy 



Us 3 



1/2 



IW| 2 dy 



1 1/2 



which proves the proposition. 



□ 



Finally we are ready to prove Theorem 14 .21 that the kernel of BS is HG(Q) © GG(Vl). 

Proof of Theorem \4~2\ First, we show that the subspace HG(Q) © GG(£l) is contained in 
the kernel of BS. Let V = V/ be in this subspace; then / is locally constant on each 
boundary component dVti and V must lie orthogonal to the boundary. 

We compute BS(Vf) using the left-translation formula © for BS, 



BS(Vf)(y) = f UVJ(x) x V^ dx 
Jn 



An 2 



L*V x f(x) dx + 2V y / L*V x f(x) ■ V y 0i dx 



Call these three terms (i), (ii), and (iii), respectively. We compute them individually, begin- 
ning with the second term. By Lemma |4~3l term (ii) becomes 



W = A / L*(V x xVJ{x))dx + -^ [ L.(V X 



f(x) x h) dx 
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Both terms on the right vanish since V/ is orthogonal to the boundary and lies in the kernel 
of curl. Hence, the second integral (ii) is zero. 

By Remark [3 .101 we may convert the gradient to x variables. 

(i) = / (L yx -i)*V x f(x) x -{L yx -i)^7 x ^Q dx 
Jn 

(Lj/x-i)* (V x /(x) x V x 4>o) dx 

The Leibniz rule for curls implies 

V/ x V0 O = V x /V0o - /(V x V0 O ) = V x /V0 O - 
After substituting this identity into (i), we apply Lemma H31 

(i) = - / (V a x /V x 0o) dx = 2 L* (fV x (j)o) dx + L* (/V^o x h) dx 
Jn Jn Jon 

Recall that / is constant on each boundary component; let fa be its value on the compo- 
nent dQi. Denote fl^ as the region "inside" <9f2j, where inside is determined opposite to the 
direction hi points. Note that fij is not necessarily a subset of Cl. Then, 

(i) = -2 / /(s)V B 0o dx + YVi / L*(Va,0oxni) dz. 
Now apply Lemma H31 again to the summed integrals: 

(i) = -2 I fVy<p dx-J^fi I ^(Vx V0 O ) dx-2j2fi \ (V" 1 )*^ dx 
Jn i Jn { i Jni 

We may now change back to gradients in terms of y variables: 

(43) (i) = -2 / fVyfodX + 2V ft I VyfodX 

Jn j Jn t 

We are now ready to compute the third integral, (iii); it will cancel the contribution of 
(i). Our first maneuver is to change gradients as above (cf. Remark [3 .101) . 

(iii) = 2V„ / {L yx -i)*V x f{x) -Vyfa dx 
Jn 

(iii) = -2V y / V«/(z) • V^i dx 
Jn 

(iii) = -2V y [ V x ■ (f(x)VM + fAfa dx 
Jn 

Recall that A<f>i = 0o + 1/87T 2 . Apply the divergence theorem to the first term, 

(iii) = -2V y I f{x)V x (f> 1 -nd{area x ) + 2V y I ffodx + 2V y / f{x)^dx 
Jan Jn Jn 



24 



JASON PARSLEY 



The last integral returns [/], the average value of / on the three-sphere, which is a constant; 
its gradient is zero. To analyze the first term, recall that / is a constant fa on each boundary 
component cftl*. Then, 

(iii) = -2V y y~]fj I V x 0i • hi d(area x ) + 2V y \ f(p Q dx 
i Jan, Jn 

(iii) = -2 V ftfy [ V x ■ V x (f)i dx + 2V y [ f((> dx, 
i JQi Jn 

where we have applied the divergence theorem. Now substitute for the Laplacian of <f>\\ 

(iii) = -2J2fiV y [ (po+strdx + 2V y [ f{x)4> dx 

Jn z Jn 

(44) (iii) = -2V/, / Vyfodx + 2 / f{x)V y <j) dx 

i Jn, Jn 

Examining equations (|43l) and (1441) . we see that (iii) is the negative of (i). Hence, 
BS(Vf) = (i) + (ii) + (iii) = 0. Therefore, the subspace HG © GG is indeed contained 
in the kernel of BS. 

Now we prove that ker BS C HG © GG. Let V G FK © HK © CG and decompose 
it as V = Vk + Vc, where V K e FK © HK and V E CG. Maxwell's fourth equation 
(Theorem 13 .31) implies that 

(45) V x BS{V K + V c ) = V K + V c -V y f O Vc • n d(area) 

J an 

If Vk 7^ 0, then it cannot possibly be cancelled by the two gradients above; in this case 

V x BS(V) 7^ and consequently BS(V) ^ 0. So it suffices to show that no curly 
gradients are in the kernel of BS. 

Assume V = V/ is a curly gradient that is in the kernel of BS. We will show then that 

V = 0. Since curly gradients are divergence-free, Maxwell's fourth equation states, for 
y e Q, that V x BS(V)(y) = V - E a , where E a is the electrostatic field d3~9i Since we 
are assuming BS(V) = on fi, the curl of BS(V) must also vanish; thus V = E a on f2. 
By Proposition I4.5L the energy of V on is no less than the energy of E a throughout the 
three-sphere. Since E a equals V inside Vt, it has no available energy left on the complement 
S 3 — f2, and so E a must be identically zero there. This implies that BS(V) = on all of 
the three-sphere, a fact which we could not conclude a priori. 

Since E a equals a gradient — Vip, the potential function ip must be locally constant on 
S 3 — £1. In particular, ip is constant on each boundary component, which implies that i*^ 
lies in HG © GG(fl). But this subspace is orthogonal to the curly gradients, hence V must 
be trivial. 

Thus we have shown that the kernel of BS cannot contain any curly gradients or knots. 
Thus it must be included in HG © GG; by the first half of the proof, the kernel is precisely 
that subspace. □ 
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Indeed, in proving Theorem |4.2[ we have shown an even stronger result, namely that no 
curly gradient can lie in the kernel of V x BS. No knot lies in this kernel due to Ampere's 
Law. Thus, the kernel of V x BS is exactly the kernel of BS. 

Theorem 4.6. The kernel ofV x BS is precisely HG{Q) © GG{Q). 

4.3. Curl of the Biot-Savart operator. An important property of the Biot-Savart operator 
is that, for certain vector fields, it acts as an inverse operator to curl. As discussed in the 
introduction, inverting curl is important to problems in plasma physics and fluid dynamics. 
Here, we precisely state when BS inverts curl. 

Proposition 4.7. (1 ) The equation V x BS(V) = V holds onVL C S 3 if and only ifV is a 
divergence-free field tangent to the boundary dQ, i.e., V G FK(Q) © HK(Q). 

(2) The equation V x BS(V) = holds onS 3 -Vl if and only ifV G FK(Q) © HK(9) © 
HG(n) © GG{Q). 

Proof. We begin with the first statement, which is far more complicated to prove than 
the second. One inclusion is immediate from Theorem 13.31 which guarantees that V x 
BS(V) = V holds onOcS 3 for a fluid knot V. 

Now for the other inclusion. Assume V x BS(V) = V holds. Then V lies in the image 
of curl, namely FK © HK © CG. The equation V x BS(V) = V holds for any fluid knot, 
and so it suffices to show that it cannot hold for a curly gradient. 

Let V = V/ G CG(tt). Theoreml33l implies that V x BS(V) = V + E a on tt, where 
E a is the electrostatic field (l39l) . Suppose on Q that V x ^^(V") = V, which implies 
E a = on n. We show that this also implies V is trivial. 

Our first step is to show that E a = outside fi. As before, express E a = —Vip. Since 
E a = on n by assumption, the potential function ij) is locally constant on Vt\ in particular 
tp is a constant tpi on each boundary component dVti. 

An exercise in Euclidean electrodynamics shows that any electrostatic field derived from 
charge on closed surface will have a jump discontinuity across the surface, but will be 
divergence-free on the interior and exterior of the surface. The same result follows from 
our definitions on the three-sphere. Applying the Leibniz rule for divergences, 



V-*pE a 



V*P-E a + ip{V-E a ) = -\E a \ 2 + i/;a, 



where a = V ■ n is defined on dQ. Now we apply the Divergence Theorem, 



s 3 -n 



E a \ 2 dx 



I 



s 3 -n 



V ■ ipE a — ipa dx 
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where hi is the normal vector on dty pointing out of S 3 — f2. Since V is a curly gradient, 
its flux through any boundary component is zero. Thus, 



(46) 



JS 3 -n a JdQi 



In Euclidean space, Gauss 's law states that the flux of an electrostatic field like E a over 
a closed surface equals the total charge enclosed by this surface. An analogous result holds 

onS 3 : 

/ E a ■ hi d(area x ) = charge Q enclosed by dVti 
Jan t 

Above when we write that charge is enclosed by dVti, we intend that the charge lies in 
S 3 — dVti. The only possible charge in the region S 3 — dVti is a surface charge g(x) = 
V(x) ■ hj which lies on other boundary components dVtj that lie outside fij. Thus, 

/ E a ■ hi d{area x ) = (±1) / cr d{area x ) = (±1) / V ■ hj d(area x ) 

The sign on the right-hand side is determined according to whether the orientations of 
and rij agree. Since V is a curly gradient, it has zero flux over any boundary component; 
thus each integral above vanishes. By equation (l46l) we conclude that 



EJ 2 dx = 0, 



s 3 -n 



and hence E a = on S* 3 — fi. 

We now have that E a = on the interiors of and S 3 — Vl. Now we apply a pillbox 
argument to a point x E dVl. Take a neighborhood of x in <9f2 and extend this into a 
"pillbox" P of small height e. Then by Gauss's law: 



E a -hd(area x ) = charge Q enc i osed = / a d(area x ). 
dP Jpudn 

But E a = on dP — dVt, so the left-hand side above is zero. Thus, for any choice of 
pillbox, J Pudn cr d(area x ) = 0. Thus a is identically zero on d£l. That implies that V is 
tangent to the boundary, so it is orthogonal to the curly gradients CG. Thus V — 0. 

Thus we have shown that the only gradients for which V x BS(V) = V holds are trivial; 
the proof of the first statement is now complete. 

The second statement is far easier to prove. For a fluid knot V, Ampere's Law guarantees 
that V x BS{V) = holds on S 3 - Q. The other subspaces HG ® GG lie in the kernel of 
VxBSby TheoremgZl Thus any V e FK®HK®HG®GG satisfies V x BS(V) = 
outside n. 

To prove the reverse implication, it suffices to show, for V a curly gradient, that V x 
BS(V) is nonzero outside f2. If it is zero, then the proof of Theorem 14.21 implies that 
V — 0. Thus the second statement holds. □ 
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4.4. Self-adjointness and image of the Biot-Savart operator. In this section we show 
that the Biot-Savart operator is self-adjoint, whether it is defined on a subdomain or on all 
of S* 3 . As a corollary, we learn something about its image. 

Proposition 4.8. The Biot-Savart operator is self-adjoint. 

Proof. Let V, W G VF(Q). We use the parallel transport version of Biot-Savart to show 

(BS(V),W) = (V,BS{W)). 



(BS(V)(y),W(y)) 



[ P yx V{x) xV y <p-W{y) dxdy 

(BS(V),W) = - [ W(y)xV y <j ) -P yx V(x)dxdy, 

by merely rearranging vectors. Since V y 4> = —P yx V x 4>, we have 

(BS(V),W) = f W(y)xP yx V x <j ) -P yx V(x)dxdy 
Jnxn 

Now parallel translate all vectors from y to x. For y ^ —x, clearly P xy P yx is the identity 
map. 

(BS(V),W) = [ P xy W(y)xP xy P yx V x <f>-P xy P yx V(x)dxdy 

(BS(V), W) = [ \ [ P xy W(y) x V X <P dy] ■ V(x) dx 
Jn Uq 

(BS(V),W) = J BS{W){x)-V{x)dx 
(BS(V),W) = (BS(W)(x),V(x)} □ 



Corollary 4.9. For Q C S 3 , the image of Biot-Savart lies in FK © HK © CG. 

For comparison, the image of BS on all of S 3 is precisely the space of fluid knots K(S 3 ). 

Proof. Since BS is self-adjoint, its image must be orthogonal to its kernel. In more de- 
tail, let W be the component of BS(V) that lies in EG © GG, the kernel of BS. Then 
BS(W) = and self-adjointness implies that (BS(V),W) = (V,BS(W)) = 0. But, 
(BS(V), W) = (W, W). SoW must be trivial, and thus the image of BS never has a 
nonzero HG © GG component. □ 
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